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Abstract 

We report about stability conditions for static, spherically symmetric objects that share the 
essential features of mass varying neutrinos in cosmological scenarios. Compact structures of 
particles with variable mass are held together preponderantly by an attractive force mediated by a 
background scalar field. Their corresponding conditions for equilibrium and stability are given in 
terms of the ratio between the total mass-energy and the spherical lump radius, M/R. We show that 
the mass varying mechanism leading to lump formation can modify the cosmological predictions 
for the cosmological neutrino mass limits. Our study comprises Tolman-Oppenheimer-Volkoff 
solutions of relativistic objects with non-uniform energy densities. The results leave open some 
questions concerning stable regular solutions that, to an external observer, very closely reproduce 
the preliminary conditions to form Schwarzschild black holes. 
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I. INTRODUCTION 



Theoretical issues concerning the understanding of physical processes where neutrinos 
are produced, scattered, or absorbed have been investigated through the years. In parallel, 
strong arguments about the neutrino contribution to the energy density of the Universe 
have been constructed on very-well understood physical basis. The role of neutrinos in the 
quintessence cosmology indeed comprises a major and still open issue that also involves the 
complete understanding of dark matter and dark energy, and the mechanism responsible for 
the onset of the accelerating phase of the Universe [l 12 1 . In particular, models of mass vary- 



ing neutrinos (MaVaN's) 13H20j coupled to a dark energy light scalar field component have 



been considered for explaining the emergence of an accelerating phase in recent cosmological 
times. 

The simplest realization of the MaVaN mechanism consists in writing down an effective 
potential which, in addition to a scalar field dependent term, contains a term related to the 
neutrino energy density. In the MaVaN scenarios the neutrinos remain essentially massless 
until recent times. When their mass eventually grows close to its present value, they form 
a non-relativistic (NR) fluid and the interaction with the scalar field stops its evolution. 
The potential energy of the dark energy component becomes the Universe's dominant con- 
tribution and the cosmological acceleration ensues. However, MaVaN models face stability 
problems for the most of choices of neutrino-scalar field couplings and scalar field poten- 
tials. Actually, when neutrinos become NR. This instability is characterized by a negative 
squared speed of sound for an effective neutrino-dark energy coupled fluid, and results in the 



exponential growth of small scale modes 



21 1 . A natural interpretation for this is that the 



Universe becomes inhomogeneous with the neutrinos forming denser structures or lumps. 
Effectively, the scalar field could mediate an attractive force between neutrinos leading to 
the formation of neutrino nuggets. 

For the scenario of structure formation, we obtain static, spherically symmetric solutions 
of the Einstein equations that describe hypothetical neutrino lumps. As a first approach, we 
study the Schwarzschild's so called interior solution for a sphere of incompressible fluid of 
constant density p and pressure p which drops from its central value to zero at the surface. 
We determine the corresponding mass defect in order to obtain the stability conditions for 
neutrino lumps. We also verify how the mass varying mechanism, as a mechanism to form 
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compact structures, can modify the cosmological predictions for neutrino masses. Due to 
the adiabatic conditions that we have set up, the thermodynamic pressure p and its explicit 
dependence on r are simply given in terms of p(r), and can be explicitly computed from 



the Tolman-Oppenheimer-Volkoff (TOV) equations for the hydrostatic equilibrium 22l. |23|. 
It justifies the assumption of an explicit analytical dependence of the variable mass on the 
radial coordinate of the massive body, m(r). Actually, at our approach the neutrino mass 
behaves like an additional extensive thermodynamic degree of freedom. 

The stability analysis reveals that our static solutions become dynamically unstable for 



modified Buchdahl limits 



24J ] set by the ratio between total mass-energy and lump radius, 



M/R. We also find regular solutions that to an outside observer resemble Schwarzschild 
black holes in formation. It remains an open question whether such stable solutions do 
exist. 

In section II, after briefly explaining the mass varying mechanism in the standard Fried- 
man - Robertson - Walker (FRW) cosmological scenario, we establish the connection among 
the energy density p, the particle variable mass m and the particle density n, for an adi- 
abatic system of mass varying particles. It results in a completely closed relation through 
which one can quantify the role of the mass varying mechanism in the study of structure 
formation 25|. We also summarize the procedure for obtaining the mass defect as a function 



of the lump radius and for performing the analysis of stability of equilibrium configurations. 
In section III we obtain the conditions for neutrino lump formation in order to verify how 
the cosmological predictions for the neutrino masses could be changed. Our analysis is ex- 
tended to TOV solutions for matter lumps with non-uniform energy densities. We draw our 
conclusions in section IV by summarizing our findings and discussing their implications. 



II. EQUILIBRIUM AND STABILITY IN THE MASS VARYING SCENARIO 



The influence of the spectrum of dynamical masses on thermodynamic of interactions are 



reported in the following. In particular, we shall investigate the effects on ste 
of neutral matter like, for instance, dark matter and sterile/ active neutrinos 



lar structures 



20 



26|]. That 



is the reason for which we shall accomplish the idea of mutual interaction among neutral 



particles driven 



mechanisms 13, 



the coupling to a scalar field through previously quoted mass varying 
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A. The mass varying mechanism 



In the usual mass varying scenario, mass varying particles are coupled to a light scalar 
field which is identified with the dark sector. Presumably, the particle mass m v has its 
origin on the vacuum expectation value (VEV) of the scalar field so that its behaviour is 
governed by the dependence of the scalar field on the scale factor. Given a particle statistical 
distribution f(q), in the flat FRW cosmological scenario, the corresponding energy density 
and pressure can be expressed by 



where q = T is the background temperature at present, and we have assumed present- 
day values with a — 1. 

Simple mathematical manipulations allow one to demonstrate that 

m W0§^ = OM) - 3p(a,0). (2) 

and, from the dependence of p on a, one can obtain the energy-momentum conservation for 
the mass varying fluid, 

p(o, </>) + 3F(p(a, 0) + p(a, <£)) = ^ - , (3) 

where H = a/a is the expansion rate of the Universe and the overdot denotes differentiation 
with respect to time (' = d/dt). The coupling between relic particles and the scalar field 
as described by Eq. (|2j) is effective when particles are NR, i. e. ~^~7^y — n ( a ) <X- a ~ A 



13Hl5|. In opposition, as long as particles are relativistic (T(a) = T /a » m(<f>(a))), the 
decoupled fluid should evolve adiabatically since the strength of the coupling is suppressed 
by the relativistic pressure increasing (p ~ 3p). The mass varying mechanism is essentially 
expressed by Eq. which translates the dependence of m on into a dynamical behaviour. 
In case of neutrinos, the most natural way to explain the smallness of the neutrino masses 
is through the seesaw mechanism, according to which, the tiny masses, m, of the usual 
left-handed neutrinos are obtained via a very massive, M, sterile right-handed neutrino for 



which the mass is driven by the VEV of a scalar field through a Yukawa coupling la, [17 ] . 



B. Cosmon lumps 



The extremely weak interaction between dark matter and/or neutrinos and a s 



owly vary- 



ing scalar field are sometimes described in terms of the cosmon dynamics |30N34| . Cosmon 
lumps are bound objects, for which the scalar and gravitational fields combine to form non- 
linear solutions of their field equations. The action of the cosmon coupled to mass varying 
particles [y or DM) can be written as 



S 



d 4 x, 



-g 



77^77 + Ig^d^d^ + V(</>) + p(<t>) 

lD7T(jr 2 



(4) 



where G is the Newton constant, R is the Ricci curvature and p(4>) is the corresponding mass 
varying particle energy density. For static, spherically symmetric solutions, one employs 
the Schwarzschild metric given by ds 2 = —B(r)dt 2 + A(r)dr 2 + r 2 (d8 2 + sin 2 (6 l )d<y9 2 ), where 
A = B^ 1 = (1 — 2GM/r)~ 1 and r is the radial coordinate. From the field equations one can 
write 

j2 



dr ; 



\r 2B 2AJ dr \ d<f) dcf> J \ dcf> 



(P " 3p) (5) 



from which the explicit dependence of on r is obtained. 

Assuming the adiabati c appro ximation at cosmological scales [21], the cosmological sta- 



tionary condition 



13 
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17] could be applied to Eq. in order to suppress its right 



hand side. Observing that the analytical dependence of the scalar field on r should disap- 
pear far from the large concentration of matter, for r values much larger than the surface 
radius R one recovers the background flatness of a FRW Universe. In this case, by adding 
the usual cosmological time dependent terms, the Eq. (jHJ) should be rewritten as the usual 
energy conservation equation for the FRW Universe, 



0. 



(6) 



16 



17|, i. e. 



(7) 



and the stationary condition should also be valid for a static approach 

dVW + dlnm(0) _ = Q 

dcf) dcf) 

The simplest scenario for this mechanism of structure formation is thus obtained from the 
extension of the cosmological stationary condition to the static configuration described by 
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Eq. (jSj). Introducing the abovementioned Schwarzschild expressions for A(r) and B(r), and 
assuming the stationary condition, Eq. fl5]) can be simplified as 



dr 2 



+ 



2-8(M/R)(r 2 /R 2 ) ~\ d0 
r - 2(M/R)(r a /R 2 )\ cb 7 
d 2 /2(r - M)\ d0 
dr2 + \ r-2M / dr 








(r < R), 
(r > 



(8) 



in 1 



2M 

™* + (j)° 1 ut ^(2M/R 3 ) 



(r < i?), 
1 



which gives 

0(r) = 
0(r) = 

where 0o?i° ni are constants to be adjusted in order to match the boundary conditions. In 
the Newtonian limit where A(r) w 5(r) ~ 1, the constants match one each other so that 
i> 1 and the above solutions reduce to 



V^2M7^ 



— arc tanh (V (2M/i? 3 



(r > i?), (9) 



LOUt 
0,1 



W = O + ri, 

r 



(10) 



which satisfies the suppression of the dependence on r, as one is away from the concentration 
of matter. 

Since <fi depends on the radial coordinate r for generical classes of curved spaces some 
ordinary mass dependencies on r can be tested. Certainly the form of m(4>) and the equation 
of state constraint may lead to quite different scenarios. Thus, we emphasize that any 
prescription for dynamical masses is model dependent, i. e. arbitrary functions for m(r) are 
equivalent to arbitrary functions for 0(r) and m{4>). 



C. Thermodynamics of "non- inter acting" mass varying particles 

To describe the connection among the extensive quantities p, m and n of an adiabatic 
system of mass varying particles in a stellar object, the thermodynamic equations can be 
summarized by 

n !hl = p + p ' 



and 



m^ = p-3p. (12) 
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The equation of state provides the dependence of pressure on energy density and specific 
entropy as p = p(p,s). Due to the adiabatic conditions, the pressure and its explicit depen- 



dence on r is simply given in terms of p(r) 25j. It can be explicitly obtained from TOV 

equations for the hydrostatic equilibrium. 

Assuming that the equation of state for a symmetrically spherical distribution of matter 

where p = p(p(r)) is determined by a univoquous dependence of p on r, we have to find 

p(m(r),n(r),r) that satisfies the system of partial differential equations given by 

dp dp dm dpdn dp 
dr dm dr dn dr dr ' 

and Eqs. (Ili p -( ll2p . Provisorily eliminating p from Eqs. ffTT]) - (fT2]) . we obtain 

d d 

4p = m— + 3n— , (14) 
dm dn 

and consequently, once setting the variable dependence a — >■ a(r), 

p(m,n) = m 1 " 3 ^* (15) 

satisfies Eq. ffl3|) . Identifying a(r) as p(r)/p(r), one finds 

p = p(m(r),n(r),r) = m l ~ s ~pn l+ ~p , (16) 

where, for simplicity, the explicit dependence of m, n, p and p on r was omitted, the constant 
k was adjusted to have p = mn'va. the NR limit, and p = n 4//s in the ultra-relativistic (UR) 
limit. Obviously, the explicit form of p(r) takes the place of an equation of state, since p(r) 
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can be determined from hydrostatic equilibrium, i. e. from TOV equations 

The above results were applied for obtaining the stability conditions of symmetrically 
spherical lumps of mass varying particles assuming them as incompressible structures of 
constant energy density p |25j | . 

D. Equilibrium and stability 

In the relativistic domain, we have 

M = M + M B = 4vr f r 2 p{r) dr, (17) 

Jo 

where 



M = 47r/ r 2 m(r)n(r) ^ ' A{r) dr, (18) 
7 



with c = 1. 

The form of Eq. ffTTl) suggests that one interprets M as the total (mass) energy inside 
a radius R, including the rest mass-energy M , the internal energy of motion W and the 
(negative) potential energy of self-gravitation U. By following the auxiliary equation given 
by Eq. (I12p . M is reduced to M (no binding energy and no motion) if p = 0. The difference 
between the total rest mass-energy and the total energy corresponds to the (positive) binding 
energy M B which keeps the stellar structure stable, 



M B = M - M = 4vr 



R 

„2 



p(r) - m{r)n{r)y/Mr) dr. (19) 



Although it is very useful to interpret M as the total energy, the internal energy of 
motion W and the (negative) potential energy of self-gravitation U are not particularly 
useful except in the Newtonian approximation where the approximation M B « U + W is 
realistic. Otherwise, the rest mass-energy M is fundamental for the stability analysis. The 
binding energy is sometimes called mass defect and it corresponds to the energy released 
during the formation of a dense star from initially rarefied matter: a typical mechanism 
for lump formation. Because of this we have M B > for a completely stable static body 
originating from diffuse matter. 

In 25] we have shown how the mass varying mechanism could modify the equilibrium 
condition based on the mass defect criterium for relativistic stellar objects. Two interesting 
effects can be achieved from such a previous analysis: one for particles with increasing 
mass inwards to the center, and the other for particles with decreasing mass. Observing 
the mass defect for compact objects composed by particles with their mass exponentially 
increasing inwards to the center, we notice that stable configurations are tremendously 
favored. Decreasing mass inwards to the center can also form stable structures up to certain 
limiting values for M/R. The equilibrium conditions are achieved just for small values of 
M/R, conditioning the existence of neutrino lumps to compact objects with arbitrarily small 



masses. The crucial point in 



25] is that coupling with the background scalar field is relevant 



in determining the limits for stability. 

In the absence of mass varying mechanisms the rest mass-energy M should be written 
as M = mN, where N is the total number of particles inside the radius R, 

N = 4tt r 2 n(r)^jA{r) dr. (20) 

J 



S 



It is important to notice that, for the static approach, the stability analysis is the same as 
when we treat mass varying particle systems with the thermodynamics governed by Eq. ffl6l) . 
As previously quoted 25| , the only difference is that the analysis has to be performed in terms 
of M (in confront with M) and not in terms of N. It is translated by the correspondence 



between M and M in terms of dM/dM . For stable equilibrium configurations [27|, |28 ] 
one always has dM/dM < 1. The stability condition is reduced to dM/dN < m for the 
case where the particle masses does not depend on the radial coordinate, m{r) = m = m . 
Because of this, the usual general relativity analysis of stability is sometimes performed 
in terms of M versus N. In case of mass varying particle systems it causes substantial 
deviations from the accurate scenario described in terms of M versus M . 

The relation between the total (mass) energy M and the rest mass-energy M and its cor- 
respondence with the precedent analysis in terms of the total number of particles N bounded 
by the surface of radius R was obtained in 25j. The stable equilibrium configurations are 
achieved when dM/dM < 1, and not when dM/dN < 1. 



III. NEUTRINO LUMPS WITH NON-UNIFORM ENERGY DENSITY 



In the usually proposed MaVaN scenarios 



3, 3, 17, 3, 26] the 



neutrinos remain essen- 



tially massless until recent times. When their mass eventually increases close to its present 
value and the interaction with the background scalar field almost ceases lfJllTJ]- The energy 
of the scalar field becomes the dominant contribution to the energy density of the Universe. 
Cosmological acceleration ensues. For the coupled neutrino-scalar field fluid the squared 
speed of sound may become negative - a signal of instability [2]]]. The natural interpretation 
of this is that the Universe becomes inhomogeneous with neutrino overdensities subject to 
nonlinear fluctuations which eventually collapses into compact neutrino lumps. 

Unlike photons and baryons, cosmological neutrinos have not been observed, so argu- 
ments about their contribution to the total energy density of the Universe are necessarily 
theoretical. The question posed here is if the fraction corresponding to the neutrino energy 
contribution is due to an isotropic and homogenous NR (p ~ 0) distribution of particles, M , 
or if it is due to a gas of weakly interacting particles (p > 0) forming lumps in an expanding 
Universe, M. The ratio given by M /M for different scenarios of mass varying mechanism 
lead to a reasonable estimative for corrections on the absolute values of the neutrino masses. 
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Current cosmological data constrain the number of active neutrino flavours as well as 
the sum of their masses to m ui < 0.75 eV at 95% c.l. [36] . This constraint does not 

i=i 

agree with the Heidelberg-Moscow bounds arising from of neutrinoless double beta decay 



which sets the limit V^m^ > 1.2 eV at 95% c.l. 37], which, however, are in agreement 
with the CMB data analysis of the WMAP results that sets 0.7 eV at 95% c. 1. for each 



neutrino species (and 2.0 eV in total) 



Actually, improvement on experimental data are 



expected to be sensitive to the effects of a finite sum of neutrino masses as small as 0.06 eV 



36. 



39], the lower limit arising by neutrino oscillation experiments that set Am 2 ~ 7.5 - 
8.7 x 10~ 5 eV 2 (2a) for solar neutrinos, and Am 2 ~ 1.7 - 2.9 x 10~ 3 eV 2 (2a) for atmospheric 
neutrinos. 

Meanwhile, all the cosmological neutrino mass predictions are performed for an isotropic 
and homogenous NR distribution of particles, which in the total absence of internal interac- 
tions and gravitational forces result in a total energy (per equivalent lump) equal to M . For 
the case where the total neutrino number N is assumed to be conserved, one can define the 
apparent (measured) value for neutrino masses as (fi) = M/N and the realistic (expected) 
corresponding value as (m) = M /N so that 

Let us then assume that the energy density of the Universe involves a gas of weakly interact- 
ing particles (neutrinos) conceived by the previously considered static, spherically symmetric 
solution of the Einstein equation with a constant energy density. Independently of the an- 
alytical characteristic of the mass varying mechanism, we assume that the neutrino mass is 
equal to m on the surface of the lump, where r = R and p = 0. 

In the Fig. [T] we represent the ratio (m) / (/j,) for different mass varying scenarios in order 
to obtain the realistic value of neutrino masses in comparison with the predicted values. It 



corresponds to the stability curves obtained in [25] . 

Obviously, the smaller is the spherical radius, the more isotropic and homogeneous is 
the cosmological neutrino distribution. For lumps in which the neutrino mass decreases 
from the surface to the center, we have a lower limit for (m)/(fi). It coincides with the 
upper limit for the radius of stable equilibrium configurations. For the other cases where 
the neutrino mass increases inwards to the center, we have completely stable equilibrium 
configurations for which (m) / (fi) decreases up to vanish on the surface of radius R. From 
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the analysis of equilibrium and stability, we have noticed that an unlimited quantity of 
(mass) energy M can be compacted as a lump when the particle mass increases inwards to 
the center of the star, constituting the preliminary conditions to form a kind of black-hole 
where the limitrophe radius correspond to the event horizon. A crude interpretation of this 
is that neutrinos forming black-holes can be practically non-massive at the their surfaces. 
Although it seems paradoxical, it can be read as a natural consequence of the mass varying 
mechanism. Actually, for these compact lumps, the most part of its mass is concentrated 
in the center. Assuming that its energy contribution to the cosmological scenario and the 
total number of particles are the same (there being or there not being lumps!), the absolute 
mass value vanishes on the surface of these objects if the Buchdahl's limit is reached. 

In general, stable scenarios lead to a shift towards smaller values for the absolute neutrino 
mass predictions. Comparing two possible scenarios for cosmological neutrinos coupling with 
a background scalar field, one with an isotropic and homogenous energy distribution, and 
the other with perturbations which results in stable neutrino lumps, the presence of neutrino 
lumps in the latter case should bring the absolute neutrino mass predictions to lower values, 
even with m < 0.07 eV, which should be problematic. 

Obtaining explicit analytical solutions of Einstein's gravitational field equations is not a 
trivial issue, on account of their complicated and nonlinear character. Even for the simple 
case of static gravitational equilibrium for a spherically symmetric matter distribution with 
uniform energy density, there are only two explicit solutions of the TOV e quat ions: the 
solution for the Einstein's Universe and the Schwarzschild solution studied in [25] • 

Our purpose is to verify how relativistic structures with non-uniform energy density 
described by TOV solutions have their equilibrium and stability conditions modified by 
the mass varying mechanism. We are focused on solutions to use in investigating fluids 
with infinite density and pressure at the center of the spherical distribution of matter. We 
consider the following two well-known TOV solutions 221 ] . 
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a. Solution A 

A(r/R, k) = 

8np(r/R, k) 
Slip(r/R, k) 



1 + 2k — k 2 



1 - (r/R)(l + 2k- tf) ( (1+2K _g )(1+2K) ) 2k1k 
2k - k 2 1 3 + 5k-2k 2 k 2 



1 + 2k - k 2 (r/RY 



K 



l + K (I + 2k - k 2 )(1 + 2k) 



(r/RY 



- (r/R)^ 



(22) 



1 + 2k-k 2 \(r/R) 2 
for which = k/(1 - 2k). 

In this case, the ratio of the total mass M to the radius R cannot reach values greater 
than 1/4 if one assumes that p(0) > 3p(0), i. e. the soft Buchdahl's limit. 
b. Solution B 



A(r/R, k) 

8irp(r/R 7 k) 
8np(r/R, k) 



{2-k 2 



K 



1 



2-k 2 (r/RY 



K 



1 - (r/RY K 



(23) 



2-k 2 {r/RY {1 + kY - {I- K) 2 (r/RY 
for which M/R = (1 - k 2 /(4 - 2k 2 ). 

In this case, the ratio of the total mass M to the radius R cannot reach values greater 
than 3/14 if one assumes that p(0) > 3p(0). 



25], in the Figs. [2] and |3] we quantify the 



By following an analogous procedure as that of 
equilibrium and stability conditions for the above TOV solutions. 

A complete similarity between the two sets of solutions can be noticed. Despite being 
difficult to observe, the gravitationally unstable configurations for Solution A, in comparison 
with Solution B, are smoothly shifted toward larger values of the ratio M/R. In comparison 
with the objects with uniform energy density, the binding energy and the stability limits, in 
terms of the ratio M/R are severally attenuated. The crucial difference between constant 
and variable energy densities concerns the form of the solutions with large values of the ratio 
M/R. The attractive interaction mediated by the scalar field in both cases, A and B, are 
not sufficiently strong to compact structures with large masses. Gravity must play a more 
relevant role for these compact objects to exist. 

For the four cases of mass dependence on r where we have neutrino masses decreasing 
inwards to the center, a striking feature is the existence of neutrino lumps with arbitrarily 
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small masses. They correspond to the lower left corner of the figures, where both, the 
mass defect and the ratio M/R vanish. Because of the small masses of this objects, the 
contribution from gravity is negligible and their existence is a consequence of the attractive 
force mediated by the scalar field. The reason is that the neutrinos are essentially massless 
inside the compact structures. For the other three cases, which includes lumps of neutrinos 
with constant mass, the stability conditions can also be achieved for massive objects with 
larger values of M/R. However, the increasing of the binding energy is limited by the soft 
Buchdahl's limit. It corresponds to an exotic equation of state at the center of the object, i. 
e. p(0) > 3p(0). Differently from which we have observed in the analysis of structures with 
uniform energy density, the conditions to form structures analogous to black-holes are not 
provided here. 

With respect to the consequences of the neutrino lump formation in how the population of 
lumps can modify the cosmological predictions for the neutrino masses, we have the scheme 
of the Fig. HJ For non-uniform densities here considered, the modifications to the mass 
predictions accomplish the stability conditions for structure formation. 

IV. CONCLUSIONS 

We have analyzed the stability conditions for static, spherically symmetric solutions of 
the Einstein equations for a system of non-baryonic matter which forms stable structures due 
to attractive forces mediated by a background scalar-field (dark energy). We have supposed 
that these astrophysical objects share the key features of the mass varying mechanism in 
cosmological scenarios. In case of dark matter scenarios, the typical size of these overdensity 
fluctuations is supposed to be large, in the range of superclusters and beyond. The natural 
interpretation of these instabilities is that the Universe becomes inhomogeneous with neu- 
trino overdensities subject to nonlinear fluctuations that eventually collapses into neutrino 
lumps. In effective terms, a scalar field mediates the attractive force among the constituting 
particles leading to the formation of stable lumps. Actually, this would turn the combined 
fluid to form compact structures which behave like cold dark matter. 

Our analysis leaves open the question whether compact structures formed by neutrinos, 
as regular and stable solutions of TOV equations, indeed exist. After reporting about 
the equilibrium and stability conditions for relativistic compact structures with uniform 
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energy density, we have extended the calculations to other two sets of TOV solutions. We 
have established the connection between the population and size of neutrino lumps and the 
modification on the standard cosmological predictions for absolute values of the neutrino 
mass. 

In addition, two distinct effects could be achieved from the study of the mass defect for 
compact objects. One for particles with increasing mass inwards to the center, and the other 
for particles with decreasing mass. For compact objects composed by particles with mass 
exponentially increasing inwards to the center, stable configurations are naturally obtained. 
In particular, compact structures with a uniform energy density has the binding energy 
increasing to infinite as the ratio M/R approximates to the Buchdahl's limit. When these 
objects shrink to such size they must inevitably keep shrinking, and eventually form a black 
hole. That is not the case of compact objects with non-uniform energy densities, where 
the conditions to form structures analogous to black-holes are not completely provided. In 
the same way, up to certain boundary values for M/R, decreasing mass inwards to the 
center can lead to stable structures. Differently from the former results, the equilibrium 
conditions are achieved just for smaller values of M/R, constraining the existence of neutrino 
lumps to compact objects with arbitrarily small masses. In fact, the coupling with the 
background scalar field is relevant in determining the bounds for stability. At the same 
time, the possibility of neutrino lumps to follow their collapse essentially due to the scalar 
mediated attractive interaction has to be more carefully elaborated. 

We have assumed that the radius and mass of the compact objects which belong to the 
family of allowed TOV solutions depend on details of the dynamical formation mechanism. 
As we have quoted before, the subtleties of the study of the neutrino mass generation 
mechanism and the analytical dependence of the scalar field on the spacetime curvature are 



not mandatory. Since the Higgs sector 20|, |26|, |4l|, |42[ and the neutrino sector are possibly 



the only ones where one can couple a new standard model (SM) singlet without upsetting 
the known phenomenology, the replacement of the explicit dependence of the neutrino mass 
on the scalar field by a direct link to the radial coordinate of the curved space is acceptable. 

To end up, we emphasize that allowing for a dynamical behaviour to a scalar field asso- 
ciated to dark energy in connection with the SM neutrinos and the electroweak interactions 
may bring important insights on the physics beyond the SM. Neutrino cosmology, in particu- 
lar, is a fascinating example where salient questions concerning SM particle phenomenology 
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can be addressed and hopefully better understood. 
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FIG. 1: Modified predictions for the absolute value of neutrino masses due to the formation of 
neutrino lumps for several ad hoc mass varying particle scenarios for compact objects with uniform 
energy density The preliminary description for the dynamical mass dependence on y = r/R, 
is given on the rectangle. For all stable scenarios the absolute neutrino mass predictions are shifted 
towards smaller values. 
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FIG. 2: The mass defect, the stability curves and the stability conditions for compact objects 

(relativistic stars) with variable density described by TOV solution A. The correspondence with 

the mass varying mechanism prescription is maintained. At this point, we restrict our analysis to 

the accurate procedure which takes into account M and M , and discards N. 
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(relativistic stars) with variable density described by TOV solution B. 
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FIG. 4: Modified predictions for the absolute value of neutrino masses due to the formation of 
neutrino lumps. We plot the ratio (m) /(//), where (m) corresponds to the absolute mass value for 
neutrinos on the surface of the compact object, i. e. r = R, and p = 0. We compare the predictions 
for type A and type B TOV solutions for relativistic stars with variable energy density. We plot 
them in the same scale {M/R G (0,0.25)) in order to verify that the modifications from one to 
another are minimal. 
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